Resonance Kondo Tunneling through a Double Quantum Dot at Finite Bias 
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It is shown that the resonance Kondo tunnehng through a double quantum dot (DQD) with even 
occupation and singlet ground state may arise at a strong bias, which compensates the energy of 
singlet /triplet excitation. Using the renormalization group technique we derive scaling equations 
and calculate the differential conductance as a function of an auxiliary dc-bias for parallel DQD 
described by 5*0(4) symmetry. We analyze the decoherence effects associated with the triplet /singlet 
relaxation in DQD and discuss the shape of differential conductance line as a function of dc-bias 
and temperature. 
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I. INTRODUCTION 

Many fascinating collective effects, which exist in 
strongly correlated electron systems (metallic compounds 
containing transition and rare-earth elements) may be 
observed also in artificial nanosize devices (quantum 
wells, quantum dots, etc). Moreover, fabricated nanoob- 
jects provide unique possibility to create such conditions 
for observation of many-particle phenomena, which by no 
means may be reached in "natural" conditions. Kondo 
effect (KE) is one of such phenomena. It was found 
theoreticalljiiiS and observed experimentallySi^ that the 
charge-spin separation in low-energy excitation spectrum 
of quantum dots under strong Coulomb blockade mani- 
fests itself as a resonance Kondo-type tunneling through 
a dot with odd electron occupation Af (one unpaired spin 
S = 1/2). This resonance tunneling through a quantum 
dot connecting two metallic reservoirs (leads) is an analog 
of resonance spin scattering in metals with magnetic im- 
purities. A Kondo-type tunneling arises under conditions 
which do not exist in conventional metallic compounds. 
The KE emerges as a dynamical phenomenon in strong 
time dependent electric fiel d^ it may arise at finite 
frequency under light illuminationii*iSii^. Even the net 
zero spin of isolated quantum dot (even A^) is not an ob- 
stacle for the resonance Kondo tunneling. In this case 
it may be observed in double quantum dots (DQD) ar- 
ranged in parallel geometrjii^, in T-shaped DQDi^ii^*i&, 
in two-level single dotsi2*i^ or induced by strong mag- 
netic fieldiSiSSiSiiSS whereas in conventional metals mag- 
netic field only suppresses the Kondo scattering. The 
latter effect was also discovered experimentallji22i242^. 

One of the most challenging options in Kondo physics 
of quantum dots is the possibility of controlling the 
Kondo effect by creating the non-equilibrium reservoir of 
fermionic excitations by means of strong bias eV 3> re- 
applied between the leads^^ {Tk is the equilibrium Kondo 
temperature which determines the energy scale of low- 
energy spin excitations in a quantum dot). However, in 
this case the decoherence effects may prevent the for- 
mation of a full scale Kondo resonance (see, e.g. dis- 



cussion m Refs. |27*28*29^. It was argued in recent dis- 
putes that the processes, associated with the finite cur- 
rent through a dot with odd J\f may destroy the coher- 
ence on an energy scale T ^ Tk and thus prevent for- 
mation of a ground state Kondo singlet, so that only the 
weak coupling Kondo regime is possible in strongly non- 
equilibrium conditions. 

In the present paper we discuss Kondo tunneling 
through DQD with even A/", whose ground state is a spin 
singlet \S). It will be shown that the Kondo tunneling 
through excited triplet state \T) arises at finite eV. In 
this case the ground state is stable against any kind of 
spin- flip processes induced by external current, the deco- 
herence effects develop only in the intermediate (virtual) 
triplet state, and the estimates of decoherence rate should 
be revisited. 

As was noticed in Ref. 15, quantum dots with even J\f 
possess the dynamical symmetry 5*0(4) of spin rotator 
in the Kondo tunneling regime, provided the low-energy 
part of excitation spectrum is formed by a singlet-triplet 
(ST) pair, and all other excitations are separated from 
the ST manifold by a gap noticeably exceeding the tun- 
neling rate 7. A DQD with even A/" in a side-bound 
(T-shape) configuration where two wells are coupled by 
the tunneling v and only one of them (say, I) is coupled 
to metallic leads {L,R) is a simplest system satisfying 
this conditioni^. Such system was realized experimen- 
tally in RefUS Novel features introduced by the dynam- 
ical symmetry in Kondo tunneling are connected with 
the fact that unlike the case of conventional SU (2) sym- 
metry of spin vector S, the 5*0(4) group possesses two 
generators S and P. The latter vector describes transi- 
tions between singlet and triplet states of spin manifold 
(this vector is an analog of Runge-Lenz vector describing 
the hidden symmetry of hydrogen atom). As was shown 
in Ref. IT^ this vector alone is responsible for Kondo 
tunneling through quantum dot with even J\f induced by 
external magnetic field. 

Another manifestation of dynamical symmetry pecu- 
liar to DQDs with even Af is revealed in this paper. It is 
shown that in the case when the ground state is singlet 
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\S) and the S/T gap S ^ Tk, a Kondo resonance chan- 
nel arises under a strong bias eV comparable with S. The 
channel opens at |eT^ — (5| < Tk, and the tunnehng is de- 
termined by the non-diagonal component Jst — {T\J\S) 
of effective exchange induced by the electron tunneling 
through DQD (see Fig. 1 (right panel)). 




II. COTUNNELING HAMILTONIAN OF 
T-SHAPED DQD 

The basic properties of symmetric DQD occupied by 
even number of electrons Af = 2n under strong Coulomb 
blockade in each well are manifested already in the sim- 
plest case n = 1, which is considered below. Such DQD 
is an artificial analog of a hydrogen molecule H2. If the 
inter-well Coulomb blockade Q is strong enough, one has 
Af — ni + Ur, ni = Ur = 1, the lowest states of DQD are 
singlet and triplet and the next levels are separated from 
ST pair by a charge transfer gap ^ Q. We assume that 
both wells are neutral at ni^r = 1- Then the effective 
inter-well exchange / responsible for the singlet-triplet 
splitting arises because of tunneling v between two wells, 
/ = v'^/Q = 6. It is convenient to write the effective spin 
Hamiltonian of isolated DQD in the form 

H,^Es\S){S\+Y,ET\Tr^){Trj\^ ^ E^X^^{1) 

V A=S,Tr) 

where X^^ — |A) (A'| is a Hubbard configuration change 
operator (see, e.g.j^), Et — Es + S, i] — ±,0 are three 
projections of S* = 1 vector. Two other terms completing 
the Anderson Hamiltonian, which describes the system 
shown in Fig.l (left panel), are 



E ^..4..^... + E E i^^^'-L.^'^ + H-c) ■ (2) 

kaa AX kaa 

The first term describes metallic electrons in the leads 
and the second one stands for tunneling between the 
leads and the DQD. Here a = L,R marks electrons in 
the left and right lead, respectively, the bias eV is ap- 
plied to the left lead, so that the chemical potentials are 
IJ-FL = IJ-FR + eV, W^"^ is the tunneling amplitude for 
the well I (left), |A) are one-electron states of DQD, which 
arises after escape of an electron with spin projection a 
from DQD in a state |A). We solve the problem in a 
Schrieffer- Wolff (SW) limili^, when the activation en- 
ergies \Ea — Ex — fiFa\ and Coulomb blockade energy 
Q are essentially larger then the tunneling rate 7, and 
charge fluctuations are completely suppressed both in the 
ground and excited state of DQD. In this limit one may 
start with the SW transformation, which projects out 
charge excitations. We confine ourselves with the bias 
eV < S <^ D, where D is the width of the electrons in 



FIG. 1: Left panel: Double quantum dot in a side-bound 
configuration. Right panel: cotunneling processes in biased 
DQD responsible for the resonance Kondo tunneling. 

the leads, so the leads are considered in the SW transfor- 
mation as two independent quasi equilibrium reservoirs 
(cf.*'^). As is shown in Ref>i^, the SW transformation 
being applied to a spin rotator results in the following 
effective spin Hamiltonian 

H^nt = Y^^iJlI'S + J!I'P) ■ Soa' + jH^X^^n^^,] (3) 

aa' 

Here^ ^aa'=X)fcfc' 'Aacr^'^k'a'a' , ''^aa'=^i^k' '^laa^'^k' a' a : 

f, 1 are the Pauli matrices and unity matrix respec- 
tively. The effective exchange constants are 

^ 2 V fiFa -Es/2^ epa' -Es/2j- 

In this approximation the small differences between sin- 
glet and triplet states are neglected. Besides, J^^, ^ I 
in real DQD. 

Two vectors S and P with spherical components 

S+ = V2 {X'° + A°-i) , S- = V2 (Xoi + A-i") , 

P+^V2 (Ai^ - A'^-i) ,P-=V2 (A'^i - A-i^) . (4) 

obey the commutation relations of 04 algebra 

[S'j,S'fe] = iejkiSi, [Pj,Pk] = iejkiSi, [Pj,Sk] = icjkiPi 

(j, fc, I are Cartesian coordinates, Cjki is a Levi-Civita ten- 
sor). These vectors are orthogonal, S • P = 0, and the 
Casimir operator is -I- P^ = 3. Thus, the singlet state 
is involved in spin scattering via the components of the 
vector P. 

We use S'J7(2)-like semi-fermionic representation for S 
operators^Si^ 

S+ = V2(/t/-i + /iVo), S- = V2(/li/o + /oVi), 

S. = flh-fUf^i, (5) 

where fl_i are creation operators for fermions with spin 
"up" and "down" respectively, whereas /o stands for 
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spinless fermion22i2^. This representation can be gen- 
eralized for 5*0(4) group by introducing another spinless 
fermion fs to take into consideration the singlet state. 
As a result, the P-operators are given by the following 
equations: 



P^ = -(/o7s + /J/o). (6) 

The Casimir operator S^ + P^ = 3 transforms to the local 
constraint 

A=±,0,s 

The final form of the spin cotunncling Hamiltonian is 

Hint — Jaa' fs fs^iaij'^k'a'cr (7) 

kk' ,aa' — L,R 

+ E/ ('^Qa''5'AA' + -^aa'^AA') ■'"o-o-'Cfcao-Cfc'a'cr'/A/A' 
kk' ,aa' AA' 

where 5'^ and P'^ {d—x,y,z) are 4x4 matrices defined by 
relations - I© and = J^^ , = J^^ and J'^^ are 
singlet, triplet and singlet-triplet coupling SW constants, 
respectively. 

The cotunneling in the ground singlet state is described 
by the first term of the Hamiltonian Q , and no spin flip 
processes accompanying the electron transfer between 
the leads emerge in this state. However, the last term in 
(UI) links the singlet ground state with the excited triplet 
and opens a Kondo channel. In equilibrium this channel 
is ineffective, because the incident electron should have 
the energy 5 to be able to initiate spin-flip processes. We 
will show in the next section that the situation changes 
radically, when strong enough external bias is applied. 

III. KONDO SINGULARITY IN TUNNELING 
THROUGH DQD AT FINITE BIAS 




FIG. 2: Leading (b,d) and next to leading (c,e) parquet dia- 
grams determining renormalization of J^{a). Solid lines de- 
note electrons in the leads. Dashed lines stand for electrons 
in the dot. 

(see, e.g.fSiSLSSiSS). These processes appear in the same 
order as Kondo co-tunneling itself, and one should use 
the non-equilibrium perturbation theory (e.g., Keldysh 
technique) to take them into account in a proper way. In 
our case these effects are expected to be weaker, because 
the nonzero spin state is involved in Kondo tunneling 
only as an intermediate virtual state arising due to S/T 
transitions induced by the second term in the Hamilto- 
nian (PI, which contains vector P. The non-equilibrium 
repopulation effects in DQD are weak as well (see next 
section, where the nonequilibrium effects are discussed in 
more details). 

Having this in mind, we describe Kondo tunneling 
through DQD at finite eV < 5 within the quasi- 
equilibrium perturbation theory in a weak coupling 
regime (cf. the quasi-equilibrium approach to descrip- 
tion of decoherence rate at large eV in Ref. . To de- 
velop the perturbative approach for T > Tk we introduce 
the temperature Green's functions (GF) for electrons in 
a dot, 5a(t) = -(T.,/A(r)/|(0)), and GF for the elec- 
trons in the left (L) and right (R) lead, GL,R{k,T) — 

— {TrCL,R(j{k, t)c]^ B.ai^' ^)) ■ Performing a Fourier trans- 
formation in imaginary time for bare GF's, we come to 
following expressions: 

C^kai^n) = {i^n — ^ka + A^L,fl)~"^, 



We deal with the case, which was not met in the pre- 
vious studies of non-equilibrium Kondo tunneling. The 
ground state of the system is singlet, and the Kondo tun- 
neling in equilibrium is quenched at T S. Thus, the 
elastic Kondo tunneling arises only provided Tk ^ S 
in accordance with the theory of two-impurity Kondo 
effect^^'^5. However, the energy necessary for spin flip 
may be donated by external electric field eV applied to 
the left lead, and in the opposite limit Tk <C S the elastic 
channel emerges at eV ~ S. The processes responsible 
for resonance Kondo cotunneling at finite bias are shown 
in Fig. 1 (left panel). 

In conventional spin S — 1/2 quantum dots the Kondo 
regime out of equilibrium is affected by spin relaxation 
and decoherence processes, which emerge at eV ^ Tk 



g°M = {iuj^-ET)-\ 77 = -1,0, 1 



g°{^n)^{terr-Esr\ (8) 

with e„ = 2TrT{n+ 1/2) and tj„ = 27TT{m+ 1/3)^^. 
The first leading and next to leading parquet diagrams 
are shown on Fig. 2. Corrections to the singlet vertex 
r(ct;, 0; a;', 0) are calculated using an analytical continua- 
tion of GF's to the real axis uj and taking into account 
the shift of the chemical potential in the left lead. Since 
the electron from the left lead tunnels into the empty 
state in the right lead separated by the energy eV, we 
have to put uj = eV, = in the final expression for 
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r(a;, 0; w', 0). Thus, unlike conventional Kondo eflFect we 
deal with the vertex at finite frequency uj similarly to 
the problem considered in Ref. ^3- We assume that the 
leads remain in equilibrium under applied bias and ne- 
glect the relaxation processes in the leads ( "hot" leads) . 
In a weak coupling regime T > Tk the leading non-Born 
contributions to the tunnel current are determined by the 
diagrams of Fig. 2 b-e. 

The effective vertex shown in Fig. 2b is given by the 
following equation 




r(2'')|', A — jST jTS 1 



'LL •-'LR 



fi^kL - eV) 

(kL + I^^L- 5 



(9) 



Changing the variable e^i^ for e^.^ — eV one finds that 



FIG. 3: Irreducible diagrams contributing to RG equations. 
Hatched boxes and circles stand for triplet-triplet and singlet- 
triplet vertices respectively. Notations for lines are the same 
as in Fig. 2 



eV) ^ Jiljli,yln{D/max{{eV - S),T}). 



Here Z? ~ is a cutoff energy determining effective 
bandwidth, j/ is a density of states on a Fermi level and 
/(e) is the Fermi function. Therefore, under condition 
\eV — (5| <C max[ey, 6] this correction does not depend on 
eV and becomes quasielastic. 

Unlike the diagram Fig. 2b, its "parquet counterpart" 
term Fig. 2c contains eV + 6 in. the argument of the 
Kondo logarithm: 



^ LR \^) - •Jll •Jlr 2^ -— 



fiekL - eV) 



(10) 



At eV S ^ T this contribution is estimated as 
ri«^(e^) - JiUFn-^n {D/{eV + S)) « rfJieV). 
Similar estimates for diagrams Fig. 2d and 2e give 
(^) - JUJIlJI^r-' In^ (i?/max{.., {eV - S),T}) 



The scaling equations for Jj^^ are as follows: 



Mr 

dlnD 



dlnD 



LL'^LRy 



dlnD 



— V 
2 



jST jTS 
^LL. + '^LR,- 



-J 



ST jTS 
LL,z'^LR,z 



(13) 



One-loop diagrams corresponding to the poor man's scal- 
ing procedure are shown in Fig. 3. To derive these equa- 
tions we collected only terms ~ {j'^)^ln"'^^{D/T) ne- 
glecting contributions containing ln[Z3/(ey)]. The anal- 
ysis of RG equations beyond the one loop approximation 
will be published elsewhere. The solution of the system 
((T^ reads as follows: 



J' 



vJ^\n{D/Ty 



J. 



ST 



jST 
•^0 



i^J^ln{D/T)' 



J' 



LR 



ST-^2 



\n{D/T) 



1 



uJ^\n{D/T)' 



(14) 



JUJIlJIr^^^^ (i?/max{c.,(ey-<5),T})x 



xln (D/max{w,ey,r}. (11) 

Then ri^^^(w) < y^l'^iuj) at eV S. 

Thus, the Kondo singularity is restored in non- 
equilibrium conditions where the electrons in the left lead 
acquire additional energy in external electric field, which 
compensates the energy loss 5 in a singlet-triplet excita- 
tion. The leading sequence of most divergent diagrams 
degenerates in this case from a parquet to a ladder series. 

Following the poor man's scaling approach, we derive 
the system of coupled renormalization group (RG) equa- 
tions for The equations for LL co-tunneling are: 



d 7^ 



dlnD 



'^'^LL -JlL: 



(12) 



Here a = L, a' = L, R. One should note that the 
Kondo temperature is determined by triplet-triplet pro- 
cesses only in spite of the fact that the ground state is 
singlet. One finds from that Tk = D exp[-l/(i' J^f )]. 
This temperature is noticeably smaller than the " equilib- 
rium" Kondo temperature T/<o, which emerges in tunnel- 
ing through triplet channel in the ground state, namely 
Tk ~ T^q/D. The reason for this difference is the reduc- 
tion of usual parquet equations for Tk to a simple ladder 
series. In this respect our case differs also from con- 
ventional Kondo effect at strong bias^i, where the non- 
equilibrium Kondo temperature T* « T^q/cV arises. In 
our model the finite bias does not enter Tk because of 
the compensation eV ^ S in spite of the fact that we 
take the argument uj = eV in the vertex ((HJ. The differ- 
ential conductance G{eV,T)/Go - jjfgp (cf. Ref. i3) 
is the universal function of two parameters T /Tk and 
eV/Tk, Go = e^f-Kh: 



G/Go - ln-2 (max[(eT/ - S),T]/Tk) 



(15) 



5 




eV/T, 



K 



20 



FIG. 4: The Kondo conductance as a function of dc-bias 
eV/Tk and T/Tk- The singlet-triplet splitting 5/Tk = 10. 



Its behavior as a function of bias and temperature is 
shown in Fig. 4. It is seen from this picture that the 
resonance tunneling "flashes" at eV ~ 5 and dies away 
out of this resonance. In this picture the decoherence ef- 
fects are not taken into account, and it stability against 
various non-equilibrium corrections should be checked. 



IV. DECOHERENCE EFFECTS 

We analyze now the decoherence rate h/rd associated 
with T/S transition relaxation induced by cotunneling. 
The calculations are performed in the same order of the 
perturbation theory as it has been done for the vertex 
renormalization (see Figs. 2 and 3). The details of the 
calculation scheme are presented in Appendix. 

To estimate the decoherence effects, one should cal- 
culate the decay of the triplet state or in other terms 
to find the imaginary part of the retarded self energy of 
triplet semi-fermion propagators at actual frequency [see 
discussion before Eq. ®], h/xd = -2ImE^(a;). The 
2nd and 3rd order diagrams determining h/rd are shown 
in Fig.5(a-d). Two leading terms given by the diagrams 
of Fig.5(a,b) describe the damping of triplet excitation 
due to its inelastic relaxation to the ground singlet state. 
These terms are calculated in Appendix (see Eqs. (|A.8|I . 
(|A.10|l 'l. One finds from these equations that the relax- 
ation rate associated with ST transition is 



ST 
d 



/D) Tnax[eV,uj,TK\ 



(16) 



It should be noted that for corrections associated with 
LL (RR) diagrams (Fig. 5a), describing co-tunneling pro- 
cesses on a left (right) lead, the use of quasi-equilibrium 
technique is fully justified when the leads themselves are 
in thermal equilibrium. We are interested in the zero 
frequency damping at resonance eV ~ S. Neglecting the 
small difference between and J^^ (seei^) , we also take 
j"^ w J"^^ — J. Thus the T S spin relaxation effect 
(|16|l does not contain logarithmic enhancement factor in 
the lowest order. It is estimated as 



(17) 



L(R) 
T T 



L(R) 
T z-*^ T 



L(R) 

a) 



R(L) 

b) 



c) 



d) 



FIG. 5: Leading diagrams (a-d) for h/rd (see text). Dashed 
line in the self-energy part stands for the singlet state of a 
two-electron configuration in the dot. 



The repolulation of triplet state as a function of ex- 
ternal bias is controlled by the occupation number for 
triplet state modified by the bias eV . The latter, in turn, 
depends on the modified exchange splitting 5* given by 
solution of the equation 



5* -6 = ReY.^i6\eV,T). 
The ReT.^ (Fig.5(a,b)) is given by 



(18) 



ReY,^^g^{ijj,eV,T) = —02 



uj In 



D 



max[aj, eV, T\ 

(19) 

where 02 ~ 1 is a numerical coefficient. As it is seen, the 
perturbative equation for ReE^ is beyond the scope of 
leading-log approximation. As a result, 5*{eV) — 5 <^ 5 
and repopulation of the triplet state is exponentially 
small. The corresponding factor in the occupation num- 
ber is 



Pt(ey) = exp(-r(el/)/T) 



(20) 



The effects of repopulation become important only at 
eV ^ 5 when \8* — 5\^ 5. In that case the quasi equilib- 
rium approach is not applicable and one should start with 
the Keldysh formalism^^'-^. This regime is definitely not 
realized in conditions considered above. 

Next 2nd order contribution is the damping of triplet 
state itself given by Eqs. (|A.12|I . HA.14() . It is seen from 
these equations that this damping is of threshold charac- 
ter: 



TT 

d 



{J/Df{uo-5)e[u~5) 



(21) 



where 9{ijj) is a Heaviside step function. These processes 
emerge only at > ^, so unlike the conventional case^i 
they are not dangerous. 

Corresponding contribution to ReE^ casts the form 



R(2) 
TTT 



{5 - Lo) In 



D 



max[((5 — cj), T] 



(22) 
where 62 ~ 1- 

Next, one have to check whether the higher or- 
der logarithmic corrections modify the estimate H17|l . 
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These corrections start with the 3-rd order terms shown 
in Fig. 5(c,d). Straightforward calculations lead to 
eV {J / D)^ \n{D / eV) correction (see the first term Eq. 
IA.23|) . This leading term like the 2-nd order term orig- 
inates from T — > S" spin relaxation processes. All other 
contributions are either of threshold character, or vanish 
at — > 0. As a result, the estimate 

h/Td ^ eViJi^/D)^[l + 0{iyJln{D/{eV))] 

holds. The topological structure (sequence of intermedi- 
ate singlet and triplet states and cotunneling processes 
in the left and right lead) in perturbative corrections for 
the triplet self energy part is different from those for the 
singlet-singlet vertex (see Appendix). Namely, the lead- 
ing (ladder) diagrams for the vertex contain maximal 
possible number of intermediate triplet states, whereas 
the higher order non-threshold log-diagrams for the self 
energy part must contain at least one intermediate sin- 
glet state. As it is seen from the Appendix (Ea. (|A.18|) - 
(|A.23|l '). the higher order contributions to /mI]T(w) are 
not universal and the coefficients in front of log have so- 
phisticated frequency dependence. As a result, the per- 
turbative series for triplet self-energy part can not be 
collected in parquet structures and remain beyond the 
leading-log approximation discussed in the Section III. 
There is no strong enhancement of the 2nd order term 
in SO{A) spin rotator model in contrast to SU{2) case 
discussed in^^. As was pointed out above, the main rea- 
son for differences in estimates of coherence rate is that 
in case of QD with odd M, the Kondo singlet develops 
in the ground state of the dot, and decoherence frustrate 
this ground state. In DQD with even Af the triplet spin 
state arises only as a virtual state in cotunneling pro- 
cesses, and our calculations demonstrate explicitely that 
decoherence effects in this case are essentially weaker. 
The 3-rd order correction to -ReS is given by 

i?eS^(^)(c.)^(^)^ln^(^) (23) 

(see Appendix) . This correction also remains beyond the 
leading-log approximation. 

Thus we conclude that the decoherence effects are not 
destructive for Kondo tunneling through T-shaped DQD, 
i.e. the Tk » h/rd is vahd provided 

5{S/Df -^Tk-^S. (24) 

This interval is wide enough because S/D <^ 1 in the 
Anderson model. 

The same calculation procedure may be repeated in 
Keldysh technique. It is seen immediately that in 
the leading-log approximation the off-diagonal terms in 
Keldysh matrix are not changed in comparison with 
equilibrium distribution functions because of the same 
threshold character of repopulation processes, so in the 
leading approximation the key diagram Fig. 5b (de- 
termining L-R current through the dot) calculated in 
Keldysh technique remains the same as (jA.8|) - (|A.10|l . 



In fact, repopulation effects result in asymmetry of the 
Kondo-peak similar to that in Ref. due to the thresh- 
old character of IvciEttt (see Appendix). This asymme- 
try becomes noticeable at eV ^ S, where our quasi equi- 
librium approach fails, but this region is beyond our in- 
terest, because the bias-induced Kondo tunneling is neg- 
ligible at large biases (see Fig. 4). 

V. CONCLUDING REMARKS 

We have shown in this paper that the tunneling 
through DQD with even M with singlet ground state 
and triplet excitation divided by the energy gap S ^ Tk 
from the singlet state exhibits a peak in differential con- 
ductance at eV « S (Fig. 4). This result is in striking 
contrast with the zero bias anomaly (ZBA) at eV w 
which arises in the opposite limit, S < Tk- In the lat- 
ter case the Kondo screening is quenched at energies less 
than S, so the ZBA has a form of a dip in the Kondo 
peak (se&i^ for detailed explanation of this effect). 

In this case strong external bias initiates the Kondo 
effect in DQD, whereas in a conventional situation (QD 
with odd Af spin 1/2 in the ground state) strong enough 
bias is destructive for Kondo tunneling. We have shown 
that the principal features of Kondo effect in this specific 
situation may be captured within a quasi equilibrium ap- 
proach. The scaling equations ((T^ . ifTHl can also be de- 
rived in Schwinger-Keldysh formalism (see Refs. 28 33) 
by applying the "poor man's scaling" approach directly 
to the dot conductance^. 

Of course, our RG approach is valid only in the weak 
coupling regime. Although in our case the limitations im- 
posed by decoherence effects are more liberal than those 
existing in conventional QD, they apparently prevent the 
full formation of the Kondo resonance. To clarify this 
point one has to use a genuine non-equilibrium approach, 
and we hope to do it in forthcoming publications. 

One should mention yet another possible experimen- 
tal realization of resonance Kondo tunneling driven by 
external electric field. Applying the alternate field V — 
VacCOs{ujt) to the parallel DQD, one takes into consider- 
ation two effects, namely (i) enhancement of Kondo con- 
ductance by tuning the amplitude of ac- voltage to satisfy 
the condition \eVac — ^ Tk and (ii) spin decoherence 
effects due to finite decoherence rate^. One can expect 
that if the decoherence rate h/r Tk, 

Gpeak/Go ^ In-' (h/TTK) (25) 

whereas in the opposite limit h/r <C Tk, 



Gpeak^G{VacCOs[iut]) (26) 

is averaged over a period of variation of ac bias. In this 
case the estimate (|15f) is also valid. 

In conclusion, we have provided the first example 
of Kondo effect, which exists only in non-equilibrium 
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conditions. It is driven by external electric field in 
tunneling through a quantum dot with even number of 
electrons, when the low-lying states are those of spin 
rotator. This is not too exotic situation because as a 
rule, a singlet ground state implies a triplet excitation. 
If the ST pair is separated by a gap from other cxcitons, 
then tuning the dc-bias in such a way that applied 
voltage compensates the energy of triplet excitation, 
one reaches the regime of Kondo peak in conductance. 
This theoretically predicted effect can be observed in dc- 
and ac-biased double quantum dots in parallel geometry. 
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APPENDIX 



We calculate perturbative corrections for E(w) by per- 
forming analytical continuation of E(iw„) into upper 
half-plane of lu. The parameter of perturbation theory 
is 1/ J <C 1 where i' denotes the density of states for con- 
duction electrons at the Fermi surface. 

The 2-nd order self energies have following structure 
(the indices T and ST in exchange vertices are temporar- 
ily omitted): 



-ki)G°{iuj2,k2)g°{iujn + iuji -f ^^2) (A.l) 
I 

over ^1,^2 in accordance with standard procedure, we 
come to following expression 



(A.2) 

I 

triplet GF. This SE has to be plugged in back to a semi- 
fermionic propagator to provide a self-consistent treat- 
ment of the problem. We denote the self-energy parts 
associated with singlet /triplet and triplet /triplet transi- 
tions as Stst and Ettt respectively. 

To prevent double occupancy of singlet/triplet states 
we take the limit Re[\s,T\ 3> T in the numerator of Eq. 
(|A.2|I . As a result, Ea. HA.2|) casts the form 



D 



D 



d^2 



tanh 



A 

2r 



tanh 



2T 



tanh 



2T 



tanh 



6- A 

2T 



(.2 - £.1 - A5,: 



/D I'D 
-D J-D 



Since all spurious states are "frozen out" we can put As = 
and Xt = S = Et — Es in denominator (in the latter 



case we perform a shift At = At ~ iirT/S) and proceed 
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with the analytical continuation icj„ uj + iO~^. Without 
loss of generality we assume lo > 0. As a result, we get 

I 



where P denotes the principal value of the integral. 

We start with discussion of self-energy parts determin- 
ing the spin relaxation due to T — > 5 transitions shown 

I 



In the opposite limit T ^ u) 

Imi:Ps^{io) ~ [J^^vfr, (A.IO) 

i?e4^1?H^(j^^.)^c.ln(^) (A.ll) 

where ln7 = C = 0.577... is the Euler constant. 

Next we turn to calculation of the triplet level damping 
due to TT relaxation processes (Fig. 5a, b). According to 
the Feynman codex, we can put iJg = at the first stage 
since the population of triplet excited state is controlled 
by finite level splitting 5. The contribution from diagram 
Fig. 5a is given by 

JmE^^^^ = /mS^^^'' ~ [JqU^ (uj — 5) 9{lo — S) 

(A.12) 

ijj — 

(A.13) 



for retarded (R) self-energies 



(A.4) 



(A.5) 



(A.6) 



(A.7) 



in Fig. 5(a,b). Assuming T <^ D and neglecting temper- 
ature corrections at low temperatures w T, we get 



vfuj (A.8) 



(A.9) 

I 

Similarly for Fig. 5b, 



ImT,^^^^ = /mSyyy^ ~ (jji/)^ (w — S) 6{uj — 5) 

(A.14) 

and, with logarithmic accuracy 



ui — 
(A.15) 

The threshold character of relaxation determined by the 

Fermi golden rule is the source of asymmetry in broad- 
ening of triplet line (see the text). Now we turn to cal- 
culation of the third order diagrams S^^^ shown in Fig. 
5c,d. 



J -D J-D 



J-D J-D W-|-^2-?l 



rfa / <^6n(6)(i - ri{^i))5{iv + 6 - a - ^) 

-D J-D 



Re^^^l'^iuj) ~ (J^u)^ r d^r r dC2n(C2)(l -n(a))P— 

J-D J-D W -1-42-41 



/m4')?(a;) ~ {J^^vf f d^2S{co + 6 - Ci) ~ (^^^^^(0))' - (J^^ 

Jo J-D Jo 



i?e4'^5(u;) ~ {J^'^u)' r dii r d^2P ~ (J^^i/)'u;ln 

Jo J-D <^ + ?2 - ?1 



FIG. 6: Fourth order leading diagrams (a-f) for triplet self-energy part. 



1^1,2,3 ki_2,3 

E(3'')(zc^„) ^ jsys Y G°(ic^i,ki)G"Hcj2,-k2)G"(ic^3,k3)^°(ic^n + lCJl+ic^2)^''(jc^„ + «CJ2 + jc^3) 

^1,2,3 ki 2,3 



Evaluation of Matsubara sums gives 



D i-D i-D 



-D J-D J-D (it^n +6 - 6 - +^2 - Cl - ^2) 



(A.16) 



D pD pD 



J-D J-D («t^n +6 - 6 - Ai)(«W„ +Ci - ^2 - A2) 



(A.17) 



r 



Let us consider first the case Ai = A2 = A5 = which 
corresponds to two singlet ferniionic lines inserted in self- 
energy part. Analytical continuation leads to following 
expression for S^^) ^ 5] (3b) j^isc) T <^ iv 



3 

JST 



1 ( D 

w in I — I — Lo 

LU 



(A.18) 



where Li2{x) is a di- logarithm function'*'^. As we already 
noticed, the first log correction to /mS appears only in 
3-rd order of the perturbation theory. Thus, 



,3 J« 



Re^TSsA^) - {J^'^^) JsT ^ 



L.2 ( 



^1 c.ln^(- 



(A.19) 



\^a[J^v) ln( -^"j 



(A.20) 



ReT,TssT{<jj) 



LO In 


D 















fofj-^j/) In 



+ ... ^ + c(w - (5) In 



UJ - 



l + dij'^v) In 



r 



(A.21) 



with coefficient a,b,c,d ~ 1. This reproduces results of 
Abrikosov-Migdal theorji^S. 



We assume now that Ai = A2 = At = S. It corresponds 
to the situation when both internal semi-fermionic GF 
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correspond to different components of the triplet. Fol- lowing the same routine as for calculation of E^^^ we find 



(3) 

TTTT 



[lo - 5) In 

r 



D 



9{uj-6) 



(A.22) 



Thus, the corrections to the relaxation rate associated 
with transitions between different components of the 
triplet have a threshold character determined by the en- 
ergy conservation. 



J 



Finally, we consider a possibility when two internal 
semi-fermionic GF correspond to different states, e.g. 
Ai — Xs — 0, whereas X2 — Xt = S. Performing the 
calculations, one finds 



(3) 

TSTT 



(5 In 



{5 - uj) In 



D 








D 




D 




e{Lo- 5)^ 


— LO 


+ 


5 In 




— w In 




5 — UJ 








J 




LU 



r 



(A.23) 



Similar expression can be derived for /TOS]^ygy(w). 

Any insertion of the triplet line in diagrams Fig.5.(a-d) 
results in additional suppression of corresponding con- 
tribution for LO < eV, which, in turn, prevents the ef- 
fective renormalization of the vertex in contrast to 
the processes shown in Fig. 3. The leading correc- 



tions in the 4-th order of perturbation theory are shown 
on Fig. 6a-f. We point out that all corrections to 

and contain an additional power of the small parameter 
S/D < 1 as ^ (5. 
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